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LONG TIME STABILITY FOR SOLUTIONS OF A /3-PLANE 

EQUATION 


TAREK M. ELGINDI, KLAUS WIDMAYER 

Abstract. We prove stability for arbitrarily long times of the zero solution for 
the so-called /3-plane equation, which describes the motion of a two-dimensional 
inviscid, ideal fluid under the influence of the Coriolis effect. The Coriolis force 
introduces a linear dispersive operator into the 2d incompressible Euler equations, 
thus making this problem amenable to an analysis from the point of view of 
nonlinear dispersive equations. The dispersive operator, L\ := exhibits 

good decay, but has numerous unfavorable properties, chief among which are its 
anisotropy and its behavior at small frequencies. 
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1. Introduction 

While the study of partial differential equations arising in fluid mechanics has 
been an active field in the past century, many important and physically relevant 
questions remain wide open from the point of view of mathematical analysis. Even 
very fundamental questions such as global-in-time well-posedness remain unresolved. 
This is particularly true in the study of the long-time behavior of systems without 
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dissipation, such as the Euler equations for incompressible flow, 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


dtu + u ■ Vu+ 'Vp = 0 inM"', 


div(M) = 0 in M"-, 

|M(a:)| —5-0 as \x\ oo. 


These model the motion of a so-called ideal, incompressible fluid occupying the 
whole space M”, where u is the (vector) velocity field of the fluid and p is the (scalar) 
internal pressure. 

The most physically relevant cases to study ((n])-([0]) are dimensions re = 2 and 
re = 3. While it isn’t apparent from the system in this form, the mathematical 
analysis of (ll.l[) - (jl.3[) in two dimensions is very different from the analysis in three 
dimensions. This is easily seen when we pass to the vorticity formulation: Define 
a; ;= V X re, the curl of re (in two dimensions w = dxU 2 — dyUi). It turns out that in 
two dimensions u satisfies a simple transport equation. 


dtu + re • Vcu = 0. 


On the other hand, the so-called Biot-Savart law allows us to recover re from the 
(scalar) vorticity w using the fact that div(re) = 0 as re = V-*-(—A)“^a;. 

The fact that oj is transported by a velocity field re which is one degree smoother 
than OJ (by the Biot-Savart law) allows one to prove that in two dimensions (jl.ip - 
p.3p is globally well-posec0 in Sobolev spaces, for s > 2. However, a priori the 
norm of the solutions may grow double exponentially in time. 


l. 1. The /3-plane model. One mechanism which helps to further stabilize the 
motion of an ideal fluid is the so-called Coriolis effect, which arises when the fluid 
is described in a rotating frame of reference, in which it is seen to experience an 
additional force (called the Coriolis force). On a rotating sphere the magnitude 
of this effect varies with the latitude and is described by the Coriolis parameter 
/. A common approximation in oceanography and geophysical fluid dynamics (see 

m, for example) is to assume that this dependency is linear. In the appropriate 
coordinate system the Coriolis parameter is then given as f = fo + fy, where the 
variable y corresponds to the latitude and /3 is a parameter (from which the model 
gets its name). In two space dimensions, this is the setting of the current paper (see 
also [20]). 

More explicitly, we are interested in studying the long-time behavior of solutions 
to the so-called f-plane equation 

dtoj -\- u-Voj = /3j^a;, 

where a; : M x ^ M, re = V-’-(—A)“^w, i?i denotes the Riesz transform in the 
first space coordinate, and /3 is a real number determined by the strength of the 
Coriolis force. For simplicity, we shall take /3 = 1. Henceforth we shall abbreviate 


(1.4) 


1 


This is not known in three dimensions, where uj is not simply transported. 
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the operator on the right hand side by 

r 

Liuj := —w. 

As described above, this system models the flow of a two-dimensional inviscid fluid 
under the influence of the Coriolis effect - we have written it here in the vorticity 
formulation. The asymmetry between the horizontal and vertical motion can be 
seen clearly in the linear term Lioj on the right hand side of the equation. In more 
technical terms, this breaks the isotropy of the original Euler equation and is a core 
feature and difficulty of the model. 

From a purely hyperbolic systems point of view, i.e. using only based energy 
estimates, one can show that solutions remain small in energy norms for a time 
span whose length is inversely proportional to the size of the initial data. There 
the operator Li does not play any role since {Liuj,lo)i ^2 = 0 and since Li commutes 
with derivatives. However, one can show that such energy norms are controlled by 
the L°° norm. Since the linear part of the equation exhibits dispersive decay in L°°, 
this controllability of higher order energy norms by the L°° norm of oj hints at a 
first possible improvement through dispersive decay. Indeed, as we already observed 
in our previous work on the related SQG equation m Remark 3.3], the smallness 
of solutions can be guaranteed for a longei0 period of time. 

The Dispersive Approach. In the present article we will improve this result and 
demonstrate smallness and decay of solutions to the /3-plane equation for arbitrarily 
long times. More precisely, we will show (see Theorem l2.ip that for any M > 0 there 
exists an e > 0 such that if the initial data are smaller than e (in suitable spaces), 
then the corresponding solution to the /3-plane equation stays small and decays until 
at least time t 6 ^ . 

We approach the question of well-posedness of the ;0-plane system for small initial 
data from a perturbative point of view, i.e. we view (II.4p as a nonlinear perturbation 
of the linear equation 


(1.5) 


dttv = LiuJ. 


Due to the properties of the operator Li this puts us in the realm of nonlinear 
dispersive differential equations: solutions to (11.51) decay in at a rate of t~^. 
This rate of decay is not integrable, so obtaining good control of solutions to the 
nonlinear equation (II.4p for long periods of time is a challenge. However, it is the 
most decay one may hope for in a 2d model, so despite being anisotropic this model 
is not further degenerate in that sense (unlike the dispersive SQG equation, which 
has a lower rate of decay of t ~2 - see [TO]!. 

However, for small frequencies Li is quite singular: a bound on third order deriva¬ 
tives in an based Besov space is required for the decay given by the linear prop¬ 
agator. For the full, nonlinear /3-plane equation one is then led to the study of 
weighted norms with derivatives in order to guarantee the decay. 


2 

|log(e)l 


More precisely, for data of size e smallness is guaranteed to hold up to a time of order 

-Q-2. 
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Further difficulties arise due to the lack of symmetries of this equation. In general, 
symmetries are a convenient way to obtain weighted estimates by commuting the 
associated generating vector fields with the equation. As noted above though, dElD 
does not exhibit a rotational symmetry. Only a scaling symmetry holds, which 
merely gives control! over a radial weight in L^. This, however, is not enough to 
control the decay of solutions. Given the form of Li, proving weighted estimates on 
the equation is therefore the key difficulty tackled in this article. 

This way we are lead to a detailed analysis of the resonance structure of the 
equation (in the spirit of the works of Germain, Masmoudi and Shatah - see for 
example m)- Here it is crucial to note that the nonlinearity u-'S/oj exhibits a null 
structure that cancels waves of parallel frequencies. However, the anisotropy of Li 
makes a precise study of the geometric structure of the resonant sets difficult or even 
infeasible. 

It is also worth noting that, compared to many other nonlinear dispersive equa¬ 
tions, this problem is highly quasi-linear in the sense that the nonlinearity has two 
orders of derivatives more than the linear operator Li. 

Related Works on Rotating Fluids. The /3-plane approximation has been used in 
numerous investigations into various types of waves in rotating fluids in two or three 
spatial dimensions. For the most part the theoretical efforts have been focused on 
viscous scenarios, such as the rotating Navier-Stokes equations. The dispersion re¬ 
lation there is anisotropic, but less singulaiB than the one we study. The viscosity 
introduces a strong regularizing effect and makes the equations semilinear. In par¬ 
ticular, Strichartz estimates can then be used to construct global solutions once 
the dispersion is sufficiently strong, and to study limiting systems as the dispersion 
tends to infinity. For the rotating Navier-Stokes equations this was done by Chemin 
et al. m, m, 0) and Koh et al. ([IB])- 

Another important achievement in the study of both the rotating Euler equa¬ 
tions and Navier-Stokes equations (in 3d) is the series of works of Babin-Mahalov- 
Nicolaenco (0, |2|, 0, |3|), which focused on studying the case of fast rotation. 
They proved that as the speed of rotation increases (depending upon the initial 
data), the solution of the rotating Navier-Stokes equations becomes globally well- 
posed and the solution of the Euler equations exists for a longer and longer time. 
This was shown in many different settings, including various spatial domains such 
as three dimensional tori ([3]) with non-resonant (0, [2]) and resonant dimensions 
(0)- Apart from the more singular dispersion relation, another key difference be¬ 
tween the works of Babin-Mahalov-Nicolaenco and the present article is that we 
are able to prove almost global existence in the setting of the whole space without 
viscosity, while their method does not give such a long time of existence and, more 
importantly, global existence without viscosity seems out of reach using those types 
of methods. In contrast, our argument does not apply to the periodic case. 


^For more details we refer the reader to Remark EH below. 

^In three dimensions it is given by ■^. However, notice that even in 3d this only gives a dispersive 
decay rate of - see e.g. [IS] . 
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Related Works in Dispersive Equations. In recent years there has been a surge of 
interest in the question of dispersive effects in fluid equations. However, just as 
is the case for the classical dispersive equations such as the (Fractional) Nonlinear 
Schrodinger Equation, in many examples (e.g. m, m) the dispersion relation is 
isotropic. This greatly helps with the resonance analysis and allows for an explicit 
understanding. In addition, rotational symmetry and scaling symmetries provide 
conservation laws for the solutions in weighted spaces. 

An example of an anisotropic equation that bears closer similarity with the /3- 
plane model is the KP-I equation, recently studied by Harrop-Grifliths, Ifrim and 
Tataru |15j . At the linear level, this 2d model produces the same decay rate as 
the /3-plane equation and does not allow for a direct resonance analysis due to the 
anisotropy. However, in contrast to the /3-plane equation, the KP-I equation does 
have enough symmetries to provide control of weighted norms and the main difficulty 
is to obtain a global decaying solution. Using the method of testing by wave packets, 
the authors were able to obtain global solutions for small, localized initial data. 

1.2. Plan of the Article. We start with the presentation of our main result. The¬ 
orem \2J\ in Section [2j This is followed by a short discussion that places it in the 
relevant context. The rest of the article is then devoted to its proof. 

We begin (Section [3]) by explaining the general bootstrapping scheme used to 
prove Theorem 12.11 It proceeds through a combination of bounds on the energy, the 
dispersion and control of some weighted norms. The individual estimates employed 
here are the subject of the remaining sections. 

Firstly, an energy inequality is established in Section [H This is a classical blow¬ 
up criterion, which shows that the growth of the energy norms is bounded by the 
exponential of the time integral of the L°° norms of the velocity and vorticity. 

We go on to discuss the dispersive decay given by the linear semigroup of the 
equation in Section [5l We explain the general dispersive result (with a proof in the 
appendix) and illustrate how we will make use of a more particular version, tailored 
to our bootstrapping argument in Section [3l 

The core of the article then is the proof of the weighted estimates in Sections [U [7] 
andm We commence by explaining in Section [6] the basics of our approach and give 
some insight into the difficulties we encounter. This is followed by an explanation of 
our strategy to prove the main weighted estimates /Propositions [7T1 17.21 and 17.3p . 

In Section [7] we embark on the proof with a few observations, which allow us to 
reduce to one main type of term, which is treated in Section [8j There we illustrate 
first the various techniques that come in handy and then go on to apply them in 
the later parts of this section. In particular, the main part of Section [8] is a detailed 
resonance analysis and relies on the observations and ideas laid out beforehand. 

In the appendix we gather the proofs of the full dispersive decay, of a corollary 
to our main result and of the global well-posedness for large data, as well as some 
computations that turn out to be useful. 

1.3. Notation. For j G Z we let Pj be the Littlewood-Paley projectors associated 

to a smooth, radial bump function —>• M with support in the shell ^ < 
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1^1 < 2 and its rescalings '■= ¥^(2 so that one has Pjg{C) = 93(2 ^C)9{C)- 

Derived from these by summation are the operators P<n , P>N etc. For example, 
P>Ar := '^ 2:1 >N gives the projection onto frequencies larger than > 0 in Fourier 
space (and can alternatively be written using a smooth Fourier multiplier (/?> 7 v(^) 
with value 1 on frequencies larger than N and supported on frequencies larger than 
N/2). 

We use standard notation for derivatives, but make use of the slight abbreviation 
to denote by a scalar derivative in ^ in any direction. Moreover, we let D := y/—A. 

For a set A we denote by XA its characteristic function. 

Whenever a parameter is carried through inequalities explicitly, we assume that 
constants implicit in the corresponding < are independent of it. 


2. Main Result 


The main result of this article is the following 


Theorem 2.1. For any M > 0 there exist cq > 0 and ko > 0 such that for all e < eo 
and k > ko, if loq satisfies 

\D^xu}o\^2 , \D^xu}o\^2 , kolz/fc < e, 

then there exist T > and a unique solution u{t) G C{[0,T]; P[^) to the initial 
value problem 


( 2 . 1 ) 


dtui + u ■ Vcj = Licj, 
< u = V^{-A)-^u, 
a;(0) = loq. 


Moreover, for t G [0, T] the solution decays in the sense that 

1 

(2.2) |w(t)|ioc + |u(t)|^oo + |Du(t)|^oo < ^ 

and the energy as well as a weight with two or three derivatives on the profile f{t) = 
e-P^ijj{t) remain bounded: 

(2.3) W{t)\H^ , |D\/(t)|^2 , \D^xf{t)\^2 < eP 

Remark 2.2 (Decay of Riesz transforms of u and w). Note that Du and uj only differ 
by a Riesz transform - since we obtain decay of u and oj through L?' based bounds 
on a weight on the profile, the proof of the theorem actually shows that any Riesz 
transforms of u and u decay at the same rate, respectively - see also Remark 15.31 
The only reason we include both Du and oo here is that they come up explicitly in 
the estimates. 


The question of stability of the zero solution was touched upon briefly in our pre¬ 
vious work m, where we showed that smallness of lj in for A; > 5 is propagated 
for at least a time T ~ e~^ |log(e)|~^. However, this result does not make use of any 
hirer structure of the equation than the blow-up criterion for the energy and does 
not give decay of the norm. 
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We point out that even for large data global solutions to the /?-plane equation 
(jl.4p can be constructed for uj in k > \ - see Theorem lB.il However, a priori, 
the norms of the solution could grow faster than double exponentially in time, 
which is completely different from the situation we have in the small data case. The 
proof of this result proceeds in analogy with that for the 2d Euler equation. We 
include it in Appendix [B] for ease of reference. 

About the Proof of Theorem \2.1[ The local well-posedness of (12.11) for A: > 1 is well- 
known. It can be seen by invoking the Sobolev embedding L°° for the 

exponential in the energy estimate (14.11) . for example. Thus the crucial part of the 
theorem is to obtain decay and existence on a “long” time interval [0,T]. 

For this we employ a bootstrapping argument, in which a decay assumption (similar 
to, but weaker than (12.21) 1 for some time interval gives the bounds (12.3p on that time 
interval. This in turn will yield the stronger decay (j2.2p . which allows us to prolong 
the solution and then repeat the argument by the continuity of w, as long as t <T. 

We make use of Fourier space methods to understand the detailed structure of 
the nonlinearity through a resonance analysis, the null structure being a crucial 
element. The present rate of dispersive decay, t~^, is not time integrable, which 
prevents us from proving estimates on a global time scale. More precisely, for the 
present quadratic nonlinearity, the bilinear estimates of type x —)• have 

a logarithmic growth when integrated in time. Together with a loss of derivatives 
in L°° in our weighted estimates, this forces us to use a lot of derivatives for the 
energy estimates. However, the constants involved therein grow exponentially with 
the amount of derivatives, so that we are not able to close our estimates on an 
exponentially long time interval. 

As outlined above in [ra the proof is spread out over the rest of this article: We 
begin by demonstrating how a bootstrapping scheme combines our estimates for the 
energy, dispersion and weights to yield Theorem 1 2. 11 Subsequently, these ingredients 
are proved; Section U] establishes energy estimates, after which we demonstrate in 
Section [5] how to obtain and use the dispersive estimate for the linear semigroup 
Finally, Sections [H [7] and [8] furnish the proof of the weighted estimates. □ 

Remark 2.3 (Symmetries). As mentioned in the introduction, apart from time trans¬ 
lation, the only continuous symmetry the /3-plane equation in (12.11) exhibits is a 
scaling symmetry. The former gives the conservation of w and u'mLf, 

\u{t)\L2 = |W(0)|^2 , \u{t)\L2 = |m(0)|^2 . 

As for the latter, one can see that if u}{t,x) solves (jl.4p . then so does cox{t,x) := 
X~^uj{X~^t, Xx) for any A / 0. Differentiating with respect to A at A = 1 gives the 
generating vector field S of this symmetry as S := —1 — tdt + x-V, and we have 

dtSio -|- u-VSuj + Su-Vlo = ^ 1 ( 50 ;). 

From this we obtain the conservation of Sio: \SLo{t)\j ^2 = \SiomL2. Moreover, 
through the commutation of S with the linear semigroup and the equation one 
can so obtain a bound on \x■'Vf {t )\]^2 (under the appropriate decay assumptions). 
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As an aside we note that from the proof of Theorem 12.11 it also follows that the 
Fourier transform of the profile also stays bounded in L°° (with two derivatives): 


Corollary 2.4. In the setting of Theorem, \2.1[ if one also has 
then 




< e, 


iciV(i,e) 


< 




The proof of this corollary builds on the techniques developed for Theorem 12.11 
and is given in detail in Appendix O 


3. Bootstrap 

We now demonstrate the bootstrap argument used to prove Theorem l2.ll deferring 
the details of the relevant estimates for the dispersion, energy and for weights on the 
profile function to the later sections of this article. The general approach here is a 
typical continuity argument for dispersive equations and has been used successfully 
in a plethora of other cases. 

We will prove Theorem 12.11 by showing that for M £ N there exist cq and ko 
(depending on M) such that for e < eo its claim holds for T > e“^, which can be 
made arbitrarily large. 

Bootstrapping Scheme. Given M £ N and ojq satisfying the assumptions in Theorem 
[Q we make the bootstrapping assumption that for 0 < t < e ^ 

1 

(3.1) |w(t)|icc + \u{t)\^oa + \Du{t)\^oa < 

From this we will deduce that the energy remains bounded, i.e. 

(3.2) \Lo{t)\jjk ^ £2, 
where /c £ N is chosen sufficiently large. 

Then we shall go on to conclude the bound for two derivatives on a weight on the 
profile f{t) = e~^^^u){t) 

(3.3) \D‘^xf{t)\^2<e^, 

from which it will follow that also three derivatives on a weight on the profile remain 
under control: 

(3.4) \D\f{t)\^,<e^. 

From this it then follows using the dispersive estimate that the following stronger 
decay estimate holds: 

1 

|w(t)|^oc + |u(t)|ioo + |T)u(t)|^oo < 

By a continuity argument this allows us to close our decay, energy and weighted 
estimates on the time interval [0,e“^]. 

The purpose of the present section is to show how these steps (13.21) . (|3.3I1 and 
(|3.4p can be derived by combining results established in later sections. 
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Energy Estimate. The blow-up criterion ()4.1|) in our energy estimate (Lemma 
l4Ti) gives 


rt A 


(3.5) 




where c(k) := Ci2^^. This we aim to bound by ea. 

Weighted Estimate. Combining Propositions 17.21 and 17.31 yields the estimate 

„f 1 

< e + ^ +P~^Nf’sf'^\D'^xf{s)\^2ds 

pt i 

+ 2^ f 1 + sN^~^ + sN~^~’^ + p~iNPsP] ds 

Jo {s) V / 

pt pt i 

+ 2 ^ / esN~‘^^ds+ / e^s-j-TT^ds, 

Jo Jo (sr 

with 0 < p <C 1 to be chosen later. Upon choosing N := 2sk and carrying out a 
time integration this reduces to the bound 

/■* i 

\D‘^xf{t)\^^<e + (^1+ |D2a;/(s)|^^ds 

-|- ktke^ ^€2 -|- es (1 -|- + e4^ 

An application of Gronwall’s inequality then gives 


(3.6) 




< 


,11/1 1 1 
e + ktke^ e 2 + e® -|- ei 


fPAp) 


using that t < e ^ with the exponent P{e, p) < es for p small enough. This we aim 
to bound by £ 2 . 

Assuming this bound on two derivatives and a weight we may use Proposition 17.II 
in direct analogy to conclude that 

pt i 

\D^xf{t)\^2<(^ + J ^ (l +p~^NPsP^\D^xf{s)\^2ds + kthh^ 

2 1/1 1 1 
+ ktke^ (e2 -l-es -|-e4 

where the extra term comes from the estimate on two weights. We thus obtain 


(3.7) 


|^V(t)L2 < 


6 “h kt ^ 6 S 6 2 -\- kt ^6^ (62 -|-6S -\- 


tPAp), 


Also this will be bounded by e 2 . 
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Dispersive Estimate. To close the bootstrap argument we need to show that from 
the bounds on the energy and the weights l (l3.2p . (j,S.3p and (13.4^ 1 we can actually 
deduce the stronger decay estimate 

1 

To this end we invoke Lemma 15.21 and deduce that for the choice N := 2s ^ as 
above we have 

\uj{t)\Loo = |e^^VLoc < \f\m+k + {\D^f\^2 + \D^xf\^2) 

with i = l + /i- 

To close we thus need to satisfy the condition 

2 2 /| 6 \ 1 1 

(3.8) A{^)tptA^^^v> <£4. 


Closing the Bootstrap. We are now in the position to show how the bootstrap 
can be closed. This is merely a matter of collecting the conditions established above 
and showing that they can indeed be satisfied. 

Given M G N we will establish that for e small enough the bootstrap (I3.ip - (I3.4I1 
closes on a time interval [0, e~^] once we choose A: G N large enough, thus concluding 
our proof of Theorem 12.11 

Firstly, to control the energy and two or three derivatives with a weight on the 
profile for 0 < f we require by (j3.5l) - (l3.7p that 


(3.9) 

(3.10) 

(3.11) 


e e 


-Mc{k)eS ^ I 


, - I - i i 

e + ke fce 2 62+68+64 


-Mt'S 


< 6 + 


6 + A ^6 fc 6862 +/ c 6 ^ 62(62 +68 +64 


-Mt^ 


< 62 . 


Secondly, to close the dispersive estimate (j3.ip we require by (j3.8jl that 
(3.12) A{fi)e pe p'e^ <£4. 


Claim: There exists a A > 0 such that for k = AM the above conditions are met. 


Proof. With this choice we have: 

(1) Condition ()3.9p is equivalent to < 1^ which in turn is simply the 

requirement that Mc{k)es < Clearly for any choice of M and k = AM 
this can be met by choosing 6 small enough (depending on M). 

(2) Conditions (I3.10p and (13.lip will be satisfied if we have 

2 AT 5 4 

A:£" —£86“^"® < £2, 

i.e. if ke& fc < 1. For k > 32M this will be satisfied once < 

1, which can be seen to hold for e small enough. 
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(3) To satisfy condition (j3.12p we need that A{n)e'^ 2 M(^,- 1 -^ + ^.^) ^ Hence up 
to the constant A{^) we require | ^ ^ One checks directly 

that this can be satisfied by choosing ^ (and thus also |) small enough. 

□ 

This completes the proof of Theorem 12.11 □ 


4. Energy Estimates 

To obtain energy estimates we differentiate the equation and get the following 
standard blow-up criterion: 

Lemma 4.1. Let uj be a solution of the ft-plane equation (II.4p on a time interval 
containing [0,T]. Then for any k & N we have the bound 

(4.1) |w(r)|^fc < |a;(0)|j:^fe exp ^ci 2^^ y \Du{t)\^oc, + \u;{t)\j^^ dt 

for some universal constant ci independent of k. 

Since our arguments will require a large number k of derivatives, here it is impor¬ 
tant to keep track of the size of the constants with respect to k. 


Proof. In this proof, all implicit constants are assumed to be uniform in k. 

For any 0 < j < /c, j € Nq, we let a £ Nq be a multiindex of order |a| = j and 
proceed by differentiating the equation by 0 “, multiply by d°‘ui and integrate over 
to obtain 


dt{d°‘w!{t),d'^uj{t))i 2 {d°‘{u ■ Vw), d'^uj)i 2 = {Lid°‘u}, d^uj) = 0. 


We note that div(u) = 0, so that in estimating {d°‘{u ■ Vw), 1,2 we only have to 
bound terms of the form \d^u-'V 8 nuj\^ 2 -: where |/3| = /-|- 1 , I 7 I = m and m + l + 1 = j. 
However, we point out that there are up to 2^ such terms. To bound one of them 
we use Holder’s inequality and Gagliardo-Nirenberg interpolation to conclude then 


that 




L2 


< 




u 


LP 


|V(9'^a;| 


L<i 


, 1 1 1 
where — I— = — 
p q 2 


<2^^\Du\^^^\Du\%^\u\% \w 


|i-i? 

'Rj 


< 2^^ \Du\j^ca l^lj^oo 

< 2 (|Zltt|^oo -|- |cj|^oo) , 


since the conditions I -- = (1 — 0)(j — 1) and m -|- 1 — - = (1 — '(?)(j — 1) imply that 

e + 'd = l = l-e + l 
Thus 

dt |w(f)|^,' < 2^^ {\Du{t)\^^ + |w(f)|^oc) \u}{t)\pjj , 


from which summation over j <k and Gronwall’s Lemma give (|4.ip . 


□ 
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5. Dispersive Estimate 

To understand the decay properties of the semigroup we adopt the point of 
view of oscillatory integrals. We write 

e^^^g{x)=[ 

Asymptotics for such integrals have been studied extensively (e.g. in [IS]). Since 
the Hessian of the phase it ^ of this integral is non-degenerate, the proof 

of the following proposition is a standard argument using stationary phase and 
Littlewood-Paley techniques and has already been hinted at in our previous work 
m - we include it in Appendix lAl for completeness’ sake. The need to control three 
derivatives of g can be seen by a scaling argument. 

Proposition 5.1. There is a constant C > 0 such that for any g G we 

hav^ 

(5-1) . 

More importantly, the goal of this section is to show how we can employ the 
weighted estimates to obtain decay of w in L°°. The main issue that arises is that 
the control of one weight in does not give integrability, i.e. \xg \]^2 does not 
control \g\]^i- Thus we cannot use (|5.ip directly. However, as will be clear from the 
details in Sections [6]ff. , it is not an easy matter to obtain other weighted estimates 
that do not grow fast, so we will work with only one weight and accept the loss this 
entails. The idea here is to avoid estimates either using energy estimates for 
high frequencies, or IT spaces with large p for low frequencies. 

Lemma 5.2. For any 0 < /i <C 1 and k we have 

(5.2) le^^VLco < |/|^3+. + 

where \ = 1 + a — -r^ and N > t) can he chosen later. 

Remark 5.3. The fact that these bounds employ based spaces shows that the 
above arguments can equally be applied to Riesz transforms of / to obtain their 
decay under the action of the semigroup In particular, this lemma shows that 
we can control w in L°° using one weight with three derivatives in (and an energy 
term). The velocity u on the other hand is one degree smoother and contains a 
Riesz transform: Recall that u{t) = e'^^*V“*“(—A)“^/(t), so this can be controlled 
using one weight with only two derivatives. 

As we will see in Sections [6ll8l our proof of Theorem 12.11 first establishes bounds 
for two derivatives on a weight of the profile /, which then gives bounds for three 
derivatives on a weight of /. In terms of the variables of the equation, we first obtain 
control of the decay of u, which then allows us to control that of uj. 

®Here Bfq denotes the homogeneous Besov space with s derivatives of the frequency localized 
pieces in U’, summed in and norm \ip\sB | (S’*^ ITiV'l_Lp)^g 2 1 , s G R, 1 < P, 9 < oo. 
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Proof of Lemma \5.^ Firstly we note that for high frequencies we may invoke the 
energy estimates via (I7.4p and obtain 

so that we may restrict to frequencies less than or equal to some > 0. 

Next we notice that by Bernstein’s inequality in two space dimensions we have 
for any // > 0 

allowing us to work with lower integrabilities than infinity. As for the dispersive 
estimate we write 

~ ~ |A*V, • PjfU, 

21 <N 2i<N 

21 <N 


as long asl + ^ = i + |. A direct computation shows that 


|e^iV-,Lp = 2'^' 






LP 


= 




LP 


Furthermore we may interpolate between ~ 1 ~ t ^ (see 

Proposition 15.ip to conclude that 




<t~^^p. 
LP'^ 


Combining this with the above estimates when r = ^ and g = 1 + /i yields 

|e^^'P<7v/Li/p 

(5-3) < t-i+lAr2/. ^ 23^2-^! 


21 <N 


1 


with - = 1 +/i - 

Now we use that by Holder’s inequality we can control by near the origin 

t2 _ _• _.. i_4.1_ 


j.'H-;vv vvc u.o^ o vv^::; y^cjo-i. ±j ±j i.j.ccn ui-ic 

and by one weight in L?' near infinity. More precisely, we have the estimates 

'JV 

.71. 


|/X|x|<H|ii+M < I/Il2 (- 7 ^ 


and 


\fX\ 


|a;|>R| — 


I — U 

i-M 

L 4/i 


l-M 

2(1+m) 


k/lL2 • 


1-M 

Optimizing R with respect to /r suggests the choice R := fj, 2 (i-i-m) ^ which leaves us 
with 

(5-4) I/Il1+m <^(/^)(l/lL2 + k/lL2), A{n) ~ n 2 (i+m)’^ . 

In particular, notice that A{fi) ^ oo as /r —^ 0. 
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Returning to the above dispersive estimate we thus see that 






Pjf\ 


Ll+M 


23 <N 


< 


(5.5) 


^-i+f^ 2^(3-|) ^ \xP,f\^,) 


23 <N 


<t 

< .-1+^ 




P<ND^—pf 




+ 


P<nD^ ’’xf 


L2 


t-^+pN^^^Aifi)NP {\P<ND^f\^, + \P<ND^xf\^,) 


Here the last inequality holds since one can combine the inequalities of Holder and 
Bernstein to show that in one has \P<nD~'^ f\i 2 ^ \P<Nf\L 2 for 0 < u < 
1 . □ 


6. Towards the Weighted Estimate 

Our goal for this and the following two sections is to provide bounds for a 
weight with two or three derivatives on the profile function, i.e. bounds for D^xf 
and D^xf in L^. This is done first for two derivatives (Propositions ES] and ES]), 
from which we can deduce estimates for the case of three derivatives (Proposition 
EH). We thus obtain first a bound for u, which then yields a bound for ui (see also 
Remark 15.3^ . 

For this we will essentially use estimates of the type x —)• L^. Informally 

speaking, in terms of our bootstrapping scheme in Section [3l every bilinear term 
then contains a decay factor of through dispersion (see (j5.ll) ) and gives back the 
weight in or an energy term, both of which are assumed to remain small. 

We start with a general discussion of our approach and give a brief overview of 
the structure of the remaining chapters and the difficulties we encounter. 


6.1. The Duhamel Term and Remarks on the Present Approach. The 

Duhamel formula for the /3-plane equation on the profile f{t) := in Fourier 

space reads 

(6.1) /(t,0 = /o(0+/ / P"'^rniC,r])f{s,C-r])f{s,r])dr]ds, 

Jo 

where the phase $ is given as 

a, = ^ ~ di _ m_ 

and the multiplier is 

\m 

or (by changing variables d i ~ d) 

i • d^ 

\i-d?' 


m{^,d) := m(^,^ - rj) 
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We note that these have a null structure which annihilates waves with parallel fre¬ 
quencies. 

As is clear from the Duhamel formula m, all estimates on /o just amount to 
assumptions on the initial data, so we may restrict ourselves to bounding the time 
integral term 

1 ^ 1 * / / e'’'''^rn{^,r])f{s,^-r])f{s,r])dr]ds 

Jo ./R2 J/2 

for / = 2, 3. 

Up to powers of |^| this leaves us with the following three main types of term^ 
to be dealt with: 

(6.2) [ f 9^ (e“'*’)m(C,r?)/(^-r7)/(r?)dr?ds, 

Jo Jr2 

(6.3) [ [ e^"‘^d^m{C,r])f{C-v)f{d)dr]ds, 

Jo Jr2 

(6.4) [ [ e^"‘^m{C,r])dJ{C-v)f{d)dr]ds. 

Jo Jr 2 

Of these, (j6.2|) is the most difficult one due to an extra factor of s from the differ¬ 
entiation (e*®^) = For (I6.3jl we may just invoke the energy estimates 

(see Section El, while (j6.4jl gives back the weight (once we notice that we may as¬ 
sume ^ and ^ — to be of comparable size, as remarked in Section El])- In particular 
we note that for these two terms the cases of two or three derivatives are completely 
analogous. It is the difficulties encountered in ()6.2I1 that guide the structure of the 
following chapters. 

More precisely, for ()6.2I) a detailed analysis of the resonance structure of the 
equation is needed. This in turn gives rise to Fourier multipliers in bilinear terms. 
However, due to the anisotropy and singularity of the linear propagator, the asymp¬ 
totics of these do not fall into the classical categories of well-studied multipliers. In 
particular, the authors are not aware of any techniques that would allow to directly 
treat the bilinear terms that so arise in a satisfactory way. 

We conclude this section with some remarks regarding the particulary difficult 
terms of the type of (j6.2l) and useful techniques for the bilinear estimates. 


6 .1.1. Difficulties in (16.21) . As we will see through the reductions in Section E this 
term is the heart of the matter and will be discussed in detail in Section [8l The 
basic issue is that we need to control the additional factor s in the time integration: 


[ [ r/)/(^ - ffififfidrjds 

Jo 4 r2 


Jo JR 2 


- v)fiv)dvds. 


6 


In what follows we may drop the explicit dependence of / on time from our notation. 
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As is well known, the typical ways of dealing with this are integrations by parts 
in either space or time (i.e. exploiting the lack of space or time resonances in various 
domains of Fourier space). In our case the anisotropy of the equation and the degen¬ 
eracy/singularity of the linear propagator make this a difficult task. In particular, 
we cannot even calculate the time resonances explicitly. 

Up to singularities when = 0, 77 = 0 or ^ = 77 we have: 

• Time resonances T := {(C)h) G x : $(^, 77 ) = 0}: They contain the 
sets {^1 = 771 = 0 } and {^1 = 771,^2 = = 1 = 772 }) for example, but their explicit 
computation is practically infeasible. 

• Space resonances S := {(^, 77 ) G x : V,j$(^,77) = 0}: These are given 
by = 277 }. 

• Spacetime resonances 7^ := SnT = {(■^, 77 ) G xM^ : 77) = 77) = 

0}: These are given by {(.^ = (0, 2 A), 77 = (0, A) : A / 0}. 

Amongst other points, through Section [7] we will see how we can reduce our study 
of this term to considering only “low” frequencies. Section [8] then explains how we 
obtain the relevant estimates through a detailed resonance analysis. More precisely. 
Sections 18.II and 18.21 lav out the available techniques for this, for which more details 
are given in Section [8.31 

6.1.2. A Remark on the Divergence Structure of the Equation. Although for most of 

what follows we will work with the above Duhamel formulation (16.ip of the /3-plane 
equation in Fourier space and use techniques inspired by the treatment of oscillatory 
integrals, it is useful to keep in mind the divergence structure of the equation: Recall 
that we have div( 77 ) = 0. In our calculations for the weighted estimate \D'^xf \^2 we 
will encounter terms of the form {u-'V{e^^^D'^xf)). It is useful to write these 

in physical space rather than Fourier space, since this way it is easier to see that a 
cancellation occurs for the highest order of derivatives on one of the profiles, i.e. 

(e“fo^ {u-V{e^^^D^xf)) ,D^xf) = {u-V{e^^*D^xf),e^^^D^xf) = 0. 

Taking into account possible localizations this shows that in such cases the null 
form does not add further derivatives. In particular, for the case of the weight with 
derivatives falling onto one of the prohles directly as in (|6.4I) . we will not lose any 
derivatives on the weight. This is important to note, for otherwise we could not 
close the estimates for the weight with a fixed amount of derivatives. 

6.1.3. A Remark on Bilinear Estimates. The basic estimate for the bilinear expres¬ 
sions that we will encounter is of type x —>• L^. Such an estimate is easily 

available if the Fourier multipliers factor onto the prohles (see Example 16.31 below). 
In the more complicated cases we will rely one of the following theorems below. 

The classical estimate by Coifman and Meyer, which gives bounds for smooth 
bilinear symbols with the appropriate decay assumptions (see IH], |9]), will be useful 
when the inputs in the Duhamel formula (j 6 .ip have comparable frequencies: 
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Theorem 6.1 (Coifman-Meyer). Let b G x M^) be smooth away from the 

origin and define the associated bilinear operator through 

Tb{f,g){fi)-= [ K^,v)fiC-v)9{v) dg. 

Assume moreover that for any multi-indices a, (3 G Nq with |a| + |/3| < 5 we have 
for some C > 0 


(6.5) 


d?d^b{C,h) <C(|e| + |ry|) 


H-l/51 


Then Tf, is a continuous bilinear map Tj, : L^xL^ ^ L'^ for 1 < p,q,r < oo satisfying 
j ^ i.e. there exists a constant K = K{p,q,r,C) such that 

( 6 . 6 ) \n{f,g)\^.<K\f\^,\g\^,. 


As noted before, the relevant setup here are estimates of the type x —>■ L^. 

On the large scale the idea is then to make use of the L°° decay properties of the 
semigroup (through the dispersive estimate (jS.ip l to control the time integrals 
in the Duhamel term () 6 .ip and its versions. 

In the case of a localization to a domain in frequency space where one of the input 
frequencies in (16.ip is much smaller than the other we shall make use of the following 
result, as presented in m (and also contained in my- 


Theorem 6.2. Assume that for 0 < s < 1 the function 6 : x —)■ M satisfies 




T OO TJS 


=: \b\j^s < oo. 


where the Littlewood-Paley projections act in the variable rj. Then the associated 
bilinear multiplier 


Tbif,g){f)-= b{f,rj)f{f,-r])g{r]) dr] 


is continuous with norm \b\j^^s as a map Tf, : x L'^ ^ L'^' for a,c satisfying 

2 


2 < a, c < 


1 - s’ 


i+i=i-£ 

a c 2 


and 1 = - + ^. 

r- ' n' 


As it turns out, the multipliers we encounter in this setup are exactly at the 
borderline s = 1 , a = oo and c = = 2 in terms of their singularities at the origin. 


Example 6.3. The obvious estimate for the bilinear part of the Duhamel term 
just amounts to the physical space estimate \u ■ Vuj\i ^2 ■ We give it here starting from 
the Fourier expression (|6.1I) to illustrate the general procedure which will later be 
employed for similar terms. We note that the null form and the phase simply factor 
on the profiles and on the integral, so that from the perspective of the Coifman- 
Meyer Theorem 16.11 we are just dealing with a bilinear multiplier with symbol 1. 
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Thus 


[ [ e“'*’m(^,r7)/(s,^-r?)/(s,r?)hr7ds 
Jo Jr2 

/■* joTi r . ri-L -- 

= / e l«l / (?-r?)e !«-''! f{s,^ -r]) ■ f{s,r])dr]ds 

Jo iK2 ?? 


i 


* • ?i 


^5- 

= / e 


^(Ve^^7)(s,e - 17) • -^(V^(l - A)-^e^i^/)(s, 77)^77 ) hs 


— I Lis . 


rt ■ 5i 


= f e^\^\^J^{u-Vuj{s))ds 

Jo 

and we conclude that 


/7 

Jo iR 2 

<-f 


e“'^7n(C, 7?)/(s, C - f?)/('S, v)dvds 
(1 


7 


7 


e hr J” ( 77 -Vw(s)) 
f {u-Vuj{s))\j ^2 ds 

Jo 

/ |7t-Va;(s)| 2^2 ds 

Jo 


by unitarity of the linear semigroup. 


In general we will follow a similar line of reasoning, the significant difference 
being the presence of more complicated multipliers. As we remarked above, because 
of the structure of the multipliers that arise we will only follow this strategy for 
frequencies that are bounded by some N > 0 (depending on time - the details being 
made explicit later). Hence one needs to check that the various multipliers have no 
singularities in the respective localizations, for one can then see that the conditions 
of the Coifman-Meyer multiplier theorem are met by the homogeneity of the terms 
involved. 


Notation. For a bilinear Fourier multiplier a{^, 77 ) we shall write 

a{c,d) ^CM icr i??71? - ??r 

to mean that satisfies the conditions of the Coifman-Meyer multiplier 

i«rkri?-»?r 

theorem EH where a,l3,'y € N. The expression |^|“ I 777 — 77 I"’' then factors onto 

the profiles and/or gives derivatives in the physical space variable x corresponding 
to 

In most cases we only explicitly state the computations required for the verification 
of the Coifman-Meyer conditions up to first order. Because of the homogeneities 
of the involved terms the further computations needed to verify the applicability of 
Coifman-Meyer’s theorem are tedious rather than difficult and are hence not always 
fully detailed. 
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7. Weighted Estimate - Reductions 


The goal for the rest of the paper is to prove the following Proposition (and its 
precursors, Propositions 17.21 and I7.3j) regarding the weighted estimate on the profile 
/ of w: 


(7.1) 




\P>ND^xf{t)\^2 ^ \P>ND^xfo\^2 + [ |w(s)|^oc |i:>^x/(s)|^2 

Jo 


+ 2 ' 


' J |a;(s)|j:^4-i-fc (|a;(s)|^oo + |u(s)|2,oo) (^1 + siV^ ^ + siV ^ ds 
rt 

, / _ _ O 7.. . , , . O 


and 


+ 2 ^ [ sN-^’^ \ujis)\l,+t ds 
Jo 


\P<ND^xf{t)\^, < \P<ND^xfo\ 

+ [ I‘^('S)Il< 

Jo 


l2 + [ N\D^xf{s)\^2\^{s)\Loods 

Jo 

\D^xf{s)\^2 ds 


(7.2) 


+ / iV (|w(s)|^oo + lu(s)|j;^oo) |w(s)|j|^i ds 
Jo 

+ J sA^ (^| m ( s)|^oo + | m ( s)|^oo |71 m ( s )|^ oo ^ |w(s)|j |^2 ds 


ft 3 

+ / p“2 |dl^ 3 :/(s)| r2 |w(s)| 2 ds 

Jo 

[ p"iA^^+^ |i:»2a;/(s)|^2 l^('S)l 2 ds 
Jo xp 


+ L p' 
+ 


where 0 < p <C 1. 


JO 

f p-^N^+P (|w(s)|^| + l^(s)l^l) l‘^(■s)lHl ds, 

^ 0 


Outline of the Proofs of Propositions \7.1\^7.^ At first we focus on proving estimates 
for two derivatives with a weight on the profile function. This is the main difficulty, 
from which the case of an additional third derivative on the weighted profile function 
follows - we explain this at the very end, in Section 18.41 

After a useful observation in Section EH for clarity we present the estimates for 
simpler terms (16.81) and (j6.4p fSections 17.21 and 17.31 respectively), which we will rely 
on later. 

Subsequently we will treat the cases of “high” and “low” frequencies separately. 
The idea is that when restricting - in a time-dependent manner - to high enough 
frequencies we can compensate for faster time growth of some of the terms through 
the boundedness of many derivatives through the Sobolev embedding - this is the 
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content of Section 17.41 and includes the proof of Proposition 17.21 In particular, no 
further manipulations will then be necessary and all multipliers directly factor onto 
the various terms, so that no special multiplier theorems are needed. 

We have thus reduced to the study of the term ()6.2p in the regime of bounded 
“low” frequencies: This is the focuse of Section [H Here we have to make use of 
the resonance and null structures in the equation to avoid too much growth in 
time. This leads us to split up Fourier space into several regions, where distinct 
resonance effects dominate. The variety of multipliers thus produced (essentially by 
integration by parts) satisfy the requirements of one of the multiplier theorems cited 
above in Remark 16.1.31 due to our localizations and our having avoided singularities 
of the multipliers involved. Due to the homogeneities of their constituent pieces, the 
conditions in the various localizations will not be too difficult to check. 

Before giving the full details in exactly the context we need them in, we lay out 
the main ideas: For low frequencies we will combine versions of a space resonance 
approach (Section 18.ip as well as time resonance methods (Section 18.21) for differ¬ 
ent regions of frequency space. The details of the splitting of Fourier space into 
appropriate subdomains and the corresponding estimates are then given in Section 

[Q 

The results are collected in Propositions 17.21 and 17.31 which give bounds for high 
and low frequencies, respectively, with two derivatives and a weight on the profile. 
Finally, Proposition 17. II contains the estimates for three derivatives and a weight on 
the profile. This concludes the proofs for the weighted estimates. □ 

7.1. A First Observation. We remark that due to the invariance of <I> under a 
change of variables 77 -fA ^ — 77 we can always arrange that I 77 I < |^ — 77 I (which 
entails |^| ~ |?~hl)' Indeed, let C ^ —>• M be a smooth cutoff function with 

supp(C) C [ 0 , 2 ] and Cl[ 0 , 1 ] = 1- Now we split the space integral in the Duhamel 
formula (16.ip by writing 

[ e“'*’777(C,r7)/(C- 77)/(77)dr7 = / e^^'^?77(C,77)(1 - C) /(C - h)/(h)^h 

4M2 V |r7| J 

+ e“‘^777(C, 77)C “ v)f{r])dr]. 

In the hrst of these terms we have I 77 I < |C — hi) hence |C| < |C “ hi ns suggested. 
In contrast, for the second term |C — hi ^ I hi) leading us to change variables to 
cr := C — h in order to get 

f e^^'^77i(C, C - cr)C ( 1 ) /(o')/(C - cr)da, 

VK-O-I/ 

since $(C,h) = ‘^(C)C “ h) = ^(C)'^)- Here we see that |(t| < |C “ i^I ns claimed. 

Conclusion: In our considerations for the weighted estimate, through 
(7.3) 

[ e“'*’? 77 (C, h)/(C - h)/(h)c^h ~ + rn) (C, h)/(C - d)f{v)X{\7^\<\i-v\}dv 

4R2 Jr2 
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we may reduce the Duhamel integral in (j6.1|) to two terms with null forms m or m 
and localization \r]\ < |^ — ?7|. 

In the rest of the article we shall make use of this localization without writing 
it explicitly every time it appears. We note that in this setting fh has a milder 
singularity than m, so often it will be enough to deal with m. 


7.2. Estimate for ()6.3jl . Since there are no derivatives on the profiles and all mul¬ 
tipliers factor (as described above), here it will suffice to invoke the energy estimate 
(j4.ip and the dispersion. Up to a Riesz transform this gives 


/ / d^rn{^,v)f{C-h)f{h)dvds 

Jo 


1 ^1^ g)f{f - r])f{ri)dri ) ds 


< 


Jo \Jm'^ 

/ \u{s)\L-°\^is)\jj2ds, 

Jo 




Ll 


and similarly for the term with m. 


7.3. Estimate for (j6.4p . The idea here is that through (j7.3p in Section 17.11 we have 
arranged that the derivative falls on the high-frequency term, which allows to 
be grouped with this derivative in order to close the estimate. The details are as 
follows. 

Writing |^| = ( |^| — \f, — vl ) +1? “ hi point out that by the divergence 

structure of the equation (Remark I6.1.2p the term with m(^, g) — r]\‘^ d^fif, — h) 
vanishes. Moreover, ~ ~ — g) • g + \g\‘^, hence 


(l^l^ - 1^ - g\‘^) m{i, g) =2 {i-g)-g 


[j-g)- g-^ 

\gf 


+ {f-g)-g-' 


< 


CM 


\^-g\"- 


When dealing with m it suffices to notice that |.^| \C-g\- 

Altogether we can thus bound two derivatives on (16.4p by 


icp r / 

Jo Jm'^ 


e^^'^m{f,g)d^f{f, - g)f{g)dgds 


f 


uj{s)\loo \f\‘^dJis,C) 




ds. 


7.4. High Frequency Estimates. For high frequencies we can make use of the 
energy estimates dnp to control the weight without resorting to any resonance 
analysis. 

Recall that we denote by P>Ar the projection onto frequencies larger than N 
in Fourier space (given by a smooth Fourier multiplier with value 1 on 

frequencies larger than N and supported on frequencies larger than N/2). The key 
point now is the following: By the Sobolev embedding we have for any /c £ N 


(7.4) 


\P>Nf\L^ < N' 


-k 


D^P>Nf 




















22 


TAREK M. ELGINDI, KLAUS WIDMAYER 


SO in particular for I derivatives 


(7.5) 


D^f 


L °° 


P>ND^f + (1 - P>N)D^f 


< 


P>ND^f 


L°° 


+ 


I 


(1 - P>N)D^f 


1/1^3+.+* +iV^ 1(1 -P>A)/li. 
<iV-" 1/1^3+.+*+iV'|/|^^. 


The moral is that by suitably splitting Fourier space we can assume the high fre¬ 
quencies to be small. 


L 2 


To estimate </?>iv(0 \C\ d^f 

ables we may assume that \C — vl 
term 


we recall from Section rrn that by changing vari- 
~ 1^1 > N, so that it is enough to bound the 




. 

Jo Jr- 


ri)P>Nf{i - v)fiv)dvds 


We then need to control the same three types of terms as in 16.11 now with an 
additional factor of ^pyN- As before, two of those are easily dealt with: When the 
weight falls on the null form m or m we proceed as in 17.21 directly invoking the 
energy estimates. In the worst case this yields 


d^rn{^,v)P>Nf{^ - r])f{rj)dr]ds 
Jo ^2 

<[ |P>ArT>^a;(s)|^2 |w(s)|^ocds. 

Jo 

Secondly, the case where the weight lands on the profile directly is dealt with as 
outlined in Section 17131 

On the other hand, term ()6.2|) (when the weight hits the phase function) requires 
some elaboration: 


The Weight on the Phase Function. As mentioned above, for high frequencies we aim 
to control the time growth from the additional factor s simply by energy estimates. 
A direct computation yields 




iC2-v2r-{Ci-vi? 

ie-r?r 

^^ 1^2 _ 2(^1 - ?? i )( g 2 - ^ 2 ) 

ler \C-vf 


so that the operator given by the symbol can be written as a sum of Riesz 

transforms in x on the one hand, and of terms with derivatives of order 2 in ^ and 
order —2 in the variable ^ — rj on the other hand. Taking into account the null form 
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as it factors onto the two profiles we thus obtain the bound 

Jo Jr2 

f s \P>nDu ;\^2 IuI^oc + S \D^ f) {e^^^D-^f)) \ 

Jo 


< 


L 2 


+ s |Z)2 (e^i*(Z)-3p>/v/) {e^^^Df))\^, ds. 


We now employ (17.51) to deduce that for the worst case 

+ |e^i*Z)-3p>/v/L2 |e"'^*^VLoc 
< \P:^nu\j^2 \Duj\]^oo + \D • 

Proposition 7.2. The high-frequency piece of the profile with weight and derivatives 

satisfies 

(7.6) 

D^xf{s )\^2 ds 

+ 2^ ^ |ca(s)|j:^ 3 +i; (|a;(s)| 2 ^oo + |'w('S)| 2 ,oo) ^1 + sN^ ^ + sN ^ ds 

+ 2^ / sN |a;(s)|j|^4-i-i: ds. 

Jo 

Proof. Note that also here we keep track of the dependence of our estimates on A; - 
all implicit constants are independent of it. 

Adding up the various pieces of the estimates from sections 171^ 17., Sl and the present 
El we obtain 


|P>/vI)V(A)| 


L 2 


< 


P>nD 


W 0 L 2 + f 

Jo 


|w(s)| 


L°° 


|P>jvZ)2a:/(t)|^2 < |P>jvD2a:/o|^2 + / |w(s)|^oc |A>^x/(s)|^ 2 ds 


(7.7) 


+ / Hs)\ |a;(s )|^2 ds 
Jo 

+ / s |P> 7 vu(s )|^2 |A>a;(s)|^oo ds 

Jo 

+ /" s|A»“^P>jva;(s)|^^ |A>^a;(s)|^2C^s- 
Jo 


Now we invoke the splitting into high and low frequencies as in (17.5p (with I = 1) 
when we have a derivative on uj in L°°, and in the same spirit we use that |P>Ar/|j /2 ^ 
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2^N This yields 

|P>ArT»^x/(i)|^2 < \P>ND^xfo\j^2 + [ |w(s)Il°° |^^^/(■s)|L2 

Jo 

+ / W{s)\Loo\uj{s)\fj2ds 

Jo 

+ 2^ y sN~’^ \uj{s)\jfk (^N~^ \uj{s)\fj4+k + N |a;(s)|j;^oo^ ds 

+ 2^ / sN~^ \uj{s)\^oo N~^ \uj{s)\fj 3 +k ds, 

Jo 

from which (17.6|) follows. 


□ 


7.5. Reduction to Low Frequencies. In the previous section we used that if ^ is 
large, then either ^ — r] or r] has to be large as well, so we can gain decay through 
Sobolev embedding. If on the other hand now ^ is small but 77 large, the same 
estimates go through. 

Conclusion: In what follows we may assume that all frequencies are bounded by 
some > 0, the time-dependence of which will be made explicit later. This means 
that we can reduce to a Duhamel formula of the type 




P<Nfo{C) + 


/7 

Jo Jr'^ 


r])P<Nf{s, ^ - r])P<Nf{s, r])dr]ds. 


For such a term we will then be able to show the following proposition: 


Proposition 7.3. The low-frequency pieee of the profile with weight and two deriva¬ 
tives satisfies 


(7.8) 


\P<ND^xf{t)\^2 < \D'^xfo\^2 + [ \D‘^xf{s)\^2 |w(s)|^oo ds 

Jo 

+ (7('S)Iloc + |u(s)|^oc) |w(s)|j;^i ds 

Jo 

+ y ■S (|w('S)licx= + |T>m(s)|^oo) \u}{s)\h 2 ds 


-T I p 2 |L)^x/(s )|^2 |w(s )|^2 ds 


+ [ P 2 A^^ (|a;(s)|^| + |u(s)|^ 2 ) |a;(s)|^i ds, 
^ 0 


where all frequeneies on the right-hand side may also he assumed to be bounded by 
N > 0 and 0 < p <C 1 (to he chosen later). 


Its proof involves the resonance analysis hinted at earlier and is carried out next 
(see Section [13]) . 
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8. Weighted Estimate - Resonance Analysis 

In this chapter we complete the proof of the weighted estimates, following the 
strategy described in the outline of the proofs of Propositions ETJES] at the beginning 
of Section^ In Sections l8.1l and [8?^ we first give a basic overview of useful techniques 
in the context of the resonance analysis for the weighted estimate. Section 18.31 
then lays out the relevant details in the case of two derivatives with a weight, thus 
concluding the proof of Proposition 17.31 Finally, in Section 18.41 we will see how one 
can obtain bounds for an additional third derivative on the weighted profile function. 
This finishes the proof of Proposition 17.11 

Many of the estimates in this section rely on the computations in Appendix [Pl 
In view of the analysis in the previous chapter, in the case of two derivatives on 
a weighted profile function we are left with estimating in the time integral 

ICP [ [ - v)f{r])dr]ds. 

Jo 

We recall that we may assume all frequencies to be bounded by some A > 0 and 
that 1771 < \^ — dl- For simplicity of notation we will not explicitly express the 
localizations, but they will be assumed implicitly in all the relevant terms that 
follow. 

8.1. The “Direct” Approach. A natural thing to do is to try to get rid of the 
extra time factor s through space resonances: Using 

( 8 . 1 ) 

IS I VrjT’l 

to integrate by parts in r] leads to terms of the type 

I ^ v^\ 

and first order derivatives in rji thereof, j, I = 1 , 2 . 

However, one computes that 

.. _ \d-m\d\" 

iv.$i \c- 2 d\\cr 

so these quotients have singularities at ^ = 0 and at ^ = 2ri and will only be useful in 
certain subdomains of Fourier space. On the other hand we point out that the null 
structure of m helps with some cancellation: since = {^ — r/)-’“-r 7 = (^ — 2 rj)-^-r], 
one order of singularity can be dealt with in some cases. However, notice that an 
integration by parts will yield a term without null structure (since the derivative 
may fall on m), in which case the singularity of ^ — 217 needs to be dealt with by 
other means. This highlights the importance of a suitable choice of region on which 
to use this approach. 

In more detail, an integration by parts in 77 using (j8.1l) leaves us with two types 
of terms: 
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(1) When the derivative in lands on one of the profiles: 

Jo JM2 \^7]^\ 

or 

f [ e**-*- m(e, 7?)^^V.cD . /(e - rj)vJ{r^)drjds. 

Jo iR2 |V,,<I>| 

This is the most delicate case to deal with, since it involves the weighted 
estimate itself - in order to be able to close our estimates we need to make 
sure that the weights on the profile come with at leasiQ two derivatives. On 
the positive side we remark that, as noted above, the null structure can help 
in controlling the singularities ^ = 2r] and = 0 in this case. 

In total, here we have a multiplier with 




ml < 


l?-2r,| 141=' 




so the singularity in ^ is only of order 2 and is thus fully compensated by 

leP- 

(2) When the derivative in lands on the Fourier multipliers: 


f [ I f{C-rj)f{r^)drjds. 

Jo ./R2 \|V,7<P| J 

Unlike the other terms, this one can be dealt with purely by energy estimates. 
For this term we cannot make use of the null structure to deal with the 
singularities arising, since it will be lost once the derivative in r] falls on it. 
Clearly we have 


V. 


|Vr,^r 


m = V 


V 




Vr,^>r 


m + Vr,<h' 






For the first term we point out that 




so that when the integration by parts produces a derivative on the phase 
quotients we have 


(8.3) 


Vr,- 




1 


2|V^$|^ 










y 

Since all frequencies are assumed to be bounded by > 0, extra derivatives just amount to a 
loss of some power of N. 
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From the computations in Appendix [D] we see that these terms are of the 
order 


1 

2 |V^<h 



I?-27?| ICI 


and 


1 









X (|17| ‘’ + I5-17I “) 


In the computations later on we will combine this with the estimates 
m\ < min{|^| , — 2r]\} \r]\~^ and |Vr,m| < |^| (see also Appendix iDl) . 


8.2. A Time Resonance Estimate. As an illustration for the estimates needed 
later on we treat here the case of pure time resonances (which will come up later on 
for Region 1 in Section 18.31) . i.e. we demonstrate how to estimate a term of the form 


Clearly 
lowing pieces: 



is$e“^m(^, 7 ])/(^ — r])f{ri)dr]ds. 


sdg so that integration by parts in time yields the fol- 



is<I>e“'^m(^, 7])/ (s, ^ — 7])f (s, r])dr]ds 

= it [ e**'^m(C, 7j)f{t, C - 7])f{t, r])d7] 

- f f e^'“^m{^,7j)f{s,^-7])f{s,7])dr]ds 

Jo Jr'^ 

- f f d)f{s,7])dT]ds 

Jo Jr'^ 

- f f se^"‘^m{^,T])f{s,^-T])dsfis,7])dT]ds. 

Jo Jr'^ 


We recall that u decays in L°° at the same rate as cu and Du (see also Remark 
15.31) . In out bootstrapping scheme all these terms will then be bounded, as the L°° 
decay gives integrability on arbitrarily long time intervals. Note that the terms with 
an extra factor of s inside the time integral are actually cubic, since by the equation 
we can substitute a bilinear term for dgf- In these cases we can take advantage of 
the decay twice. We proceed to give the relevant details. 


The first term. We note that 


r 


-1 


m{^, rj)f{t, ^ - 7])f{t, ri)drj = u-Vuj{t) 
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SO this is easily controlled: 








|t'u-Vw(t)|^2 <t\u{t)\^oo |Va;(t)|^2 • 


The second term. This is of the same order as the first term, since 


r?)/(s, ^ - r?)/(s, r])dr] 


ds= \u-'S/oj{s)\j ^2 ds. 

Jo 


The third and fourth terms. The third and fourth terms can be estimated essentially 
in the same way, so we only give the details for the third one. Prom the definition 
of the profile f{t) := e~^^^uj{t) it follows that 

dtf{t) = e~^^\-Liu;{t) + dtcoit)) = e~^^\u-Vuj{t)). 


Hence the standard estimate gives 




f [ se"'''^m{^,r])dsf{s,^-r])f{s,r])d'n 
Jo Jr^ 

= f s\u-V [e^^^{u-Vu ))\^2 ds 
Jo 

< / s\u{s)\];^oo\V {u-Vlo )\]^2 ds 
Jo 

< / s|u(s)|^oc (|u(s)|^oo |w(s)|jy2 + |T>u(s)|^co |u;(s)|^i)(is. 

Jo 


8.3. The Splitting of Fourier Space for Low Frequencies. We will use the 
following divisions of Fourier space into regions where various techniques will be 
employed to control (16.21) . the most delicate term that comes up when we estimate 
a weight on the profile. We recall fSection 17.51) that for the estimates themselves we 
may assume that all frequencies are bounded by some > 0 (which of course is not 
relevant for the splitting, but rather for the potential derivative loss), as well as the 
restriction that It;! < \ f — r]\, as discussed in Section ITTl 
We point out our use of the following 

Notation. In this section, all integrals in jy have a localization (which we may not 
write explicitly) given by the particular subsection they are found in. 


Region 1: 


(8.4) 

and 

(8.5) 


|^<l^l<100|ry| 


1^1 


10000 


< |^-?7| < 10000 I 


Region 2: 
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Region 3: 


(8.7) 


ICI > 100|r/| 


Let us now discuss how to obtain the relevant bounds for (j6.2p . As outlined further 
above, this will involve a study of the resonances (by which the above splitting of 
Fourier space is inspired) in the various domains and subsequent use of the techniques 
discussed in Sections ED and 18.21 

The result of this process is the estimate ()7.8p from Proposition 17.31 (0 < p <C 1): 


\D'^xf{t )\^2 < |T>\/o |^2 + / \D^xf{s )\^2 |w(s)Iloc ds 

Jo 

/ (|w(s)Iloc + |u(s)|^oo) |w(s)|^i ds 
Jo 

j s (|u(s)|^oo + |u(s)|ioo |Du(s)|^oc) |w(s)|j :^2 ds 

ft 3 

+ / p~ 2NP \D‘^xf{s)\^ 2 \^{s)\ ids 
Jo 

+ f p-^NP (^\ujis)\j+\u{s)\^g'^\uj{s)\fjids. 

J 0 


+ 


+ 


Here the first three integrals come from the estimates in Regions 1 and 2 via the 
Coifman-Meyer Theorem 16.11 whereas the last two integral terms are the product 
of an estimate using Theorem 16.21 in Region 3. The basic techniques for the proof 
of this estimate have been introduced above, so now we complete the steps outlined 
there by discussing the bounds for the various regions of our splitting. 


8.3.1. Region 1 - (|8.4p and (j8.5p : |^| rv-' I fj I I'Nj 1^ — ? 7 |. Here all frequencies are of 
comparable size and hence the multipliers can be controlled via Coifman-Meyer 
type estimates. The main difficulty is to deal with the space resonance ^ = 2rj. 
We thus split into two cases depending on the size of — 2r]\: Away from the space 
resonance a direct integration by parts goes through, whereas for “small” — 2p| we 
further (anisotropically) split into two subcases and use time resonances or a “good” 
direction in which to integrate by parts. These steps rely on the null structure of 
the bilinearity. 


Case 1: 

( 8 . 8 ) 


ie-2pi > 


Idl 

1000 


In this case all singularities are avoided and we will only need to integrate by parts 
in T], Indeed, recall from (j8.2|) in Section 18.II that 


Vg^l _ \il-2C\\iif 

|^_27?||er 
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Under ()8.4p . (j8.5p and ()8.8p we see directly thalU 

|V^$| 


5.CM 1- 


|V^$| 

Hence this term (of type ([T]) as described in Section [8TT|) gives back the weight D^xf 
in and ui in L°°, i.e. for it we have the bound 

1 


[ \CfdJ{s) 

Jo 


Similarly, from (18.31) and the observations in Section 18.11 we also deduce that 






^CM 


\ri\-\ 


so that the multipliers for terms of type ([2|) are of order \r]\ ^ (in the sense of 
Theorem 16.ip and we can bound them by 

rt 


Case 2: 

(8.9) 

Here we observe that 

( 8 . 10 ) 


Jo 

\i-24<2L 

'' “ 1000 

^ l^il “ T(lx) l^il 


|4>l > 




Indeed, this is gotten simply by noting that 


^ = 6 


1 


1 


1^1^ \i-ri? 


+ m 


\vf 


and invoking (18.91) and (18.41) and (18.51) . 


Proof of (j8.10p . We begin by observing that — 2r]\ < 4^ implies that \r]\ > 


14 > -mm 1^1 and thus \r]\ > \f, - r]\ > 


999 


1^1- 


Now one has 


\^-vf 




1 


which yields 


1 

\(-ri\ \ri\ 

1 






and thus 
1 


\v\\f-v\ \v\\^-r]\ 1000|C-r/|’ 


hr 


< 


< 


lOOOhll^-r/l 9991 


m 


^Recall the notation <cm for boundedness in the sense of a Coifman-Meyer multiplier from page 
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On the other hand, 


I 

1 


1 

I 

I 1 

1 


1 

UT’ 

1 

Uy> 

1 1 

1 

id^ 

le-d^ 


Id 

ic-f?i 

Id \^-v\ 


Id 1C-6 

Id IC-d 


> 


> 


1999 I I I I 

1000 I 'I 1000 I "I 
6 1 


1000 1 1000 1 

+ 


2001 |?7| 1001 |?7| 


ioir?r 

Altogether we have thus established (IS.lOp . 


□ 


This leads us to two subcases. 
Subcase A : 

( 8 . 11 ) 


leii > 


\m\ 

100 


In this case, from (IS.lOp we deduce that 


l$l > 


2|r?| 


2 • 


But then 


\m\ 


1 16%-?2r/i| ^ 2|r/|^ 161 (hi + Id) 




\v\ 


\vf " 161 

by virtue of (j8.4[) . (j8.5p and (I8.11h . It follows that ^ dcM hi and hence we can 
integrate by parts in time to obtain estimates similar to the ones before (where the 
shorthand XA denotes the localizations of this subcase): 


f [ ^66^ 

Jo ./R2 


9^4>e“‘^m(6 v)XA fiC - v)f{v)dr]ds 


( 8 . 12 ) 


/7 

Jo Jr'^ 

uf 

Jm 


(e ) "^-(6 V)XA /(? - v)f{v)dvds 


m 


e* Id d^^-^XAf{^-v)fiv)dv 


m 


e“ Id d^^—XAfi^-v)fiv)dvds 



10 iR2 

. 

Jo Jr"^ 


se“'*’ Id^ d^^^XA dsfi^ - rj)f{r])dr]ds 


m 


se“ Id d^<^—XAfi( - r])dsf{v)dvds. 

These terms can now essentially be estimated as outlined in Section [8.21 Inote that 
the additional multipliers hP are of order one, since |d ~ h ~ 3?| in Region 1). 
In particular, no weights are required and there is no derivative loss in 
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Subcase B: 


(8.13) 



In this case, integration by parts in time will not avail us. However, due to the 
constraints on and the other variables we may integrate by parts in 772 to obtain 
something useful: 


Lemma 8.1. Assume that — 2ri\ < and |^i| < Then 



(8.14) 


and 



(8.15) 


These estimates allow us to integrate by parts in 772 in this region. One notes that 
as necessary we end up with two derivatives on the weight and lose one derivative in 
the energy estimate - the arguments are completely analogous to our prior treatment 
of Case 1 for Region 1 (this is (|8.4p and (|8.5p with (|8.8p '): the integration by parts 
in space proceeds in close analogy to the corresponding observations on page [29] and 
we end up with the bound 



in this subcase. 

Proof of Lemma, \8.1\ We begin by recalling that 


(8.16) 



Note also that 



which implies 



In direct analogy we have 



and 




and thus 














LONG TIME STABILITY FOR SOLUTIONS OF A /3-PLANE EQUATION 


33 


To estimate dri 2 ^ we use that for the terms in (I8.16|) we have 1 171 ^ 2 ! > |hi| Ivl 

4 

100 


|6(% - 6)1 < w 1^2 - 61 < tm \Vi\\v\, so that 


I??i6 + 6??2-66I > 


In addition, 

Vi'n2 - 

which gives the claimed bound 

15772*^1 ^ 2 — 


\v\^ 


1 


< ^ i-fii 111 . 


\vi\\v\\m\\v\ 


V 4 100; 4 |e- 7 ?r 

For the null form we notice that with the present restrictions on ^,r] we have 


? • V 




\m\\v\ ■ 


On the one hand |6^‘'‘| = |-6% + 6^i| > since |6%l < 1^1 

|6ml > \v\\m\- On the other, |6 h~^| ^ | 6 ?? 2 | + |6??i| < l?/i| (tto + l^al) <4|r7 i||?7 |. 
To control the remaining term (I8.15jl for the integration by parts we write 


5, 


r /2 


m 

5^2 4> 


19772^1 


{dr^2'i^ 9772 4* - m 9^2^) • 


We recall from Appendix [Pl that 


9^ d> = 


’ \^-vf \C-vX 


+ - 


Vi 


r]l 


6 ®I„|4’ 


\vX \v\ 


so that 


Furthermore, 


|94d>| 


< M 


hr 

1 

I? 

3 


|9^2"i(6??)I = 7-4 -6 hP - 2??2(6 


-Tr(^!’)!" + 2 1'fll’iill'fl 




Putting these together we obtain 

m 


9 


m 


9772 $ 


< 


hr hirhi 


hrhil'hr 


,4 < 


Due to the present localization to frequencies of comparable sizes and the homo¬ 
geneities of the terms involved this yields (|8.15p . □ 
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8.3.2. Region 2 - (18.6p .' |^| <C |r/| ~ |^ —f/l- Here |^| < so that |^ — 2 t/| > \ri\ 
and we see that direct integration by parts in t] goes through, in analogy to Case 1 
of Region 1. More precisely, from 

(8-17) 1^1^ 

we get back the weight with two derivatives and lose one derivative in the energy 
estimates: once we notice that in this Region we have I < |ry| ^ 1^1 it follows 
(see also the discussion for terms of type (j2|) in Section 18.ip 

mV^ • <CM ■ 

8.3.3. Region 3 - (18.71) .' \ri\ ^ |^| ~ l'^ —??|- We have \^ — 2rj\ > |^| > \r]\, so we 

can proceed with the same integration by parts as in Region 2, given by (|8.6p . 
One can check from the localization (into “high” and “low” frequencies) that the 
multipliers arising here are not of Coifman-Meyer type. We will instead use Theorem 
16.21 to control them. By homogeneity of the constituent terms (and the frequencies 
being bounded by some > 0) one sees directly that any such multiplier b{^, rj) 
is bounded in and thus (see Lemma 5.1 in m) one also has < oo 

for 0 < s < 1. More precisely, for 0<p<Clwelets = l — p and have the 
bound \b\j^i-p < p~ 2 N^. By Theorem 16.21 the associated bilinear operator then is 
bounde41 as a map L'^ x L'p ^ L^. 

Remark 8.2. We point out that it is this region that requires us to work with two 
derivatives on the weight, rather than three. In all other regions three derivatives 
would give back three derivatives on the weight (and possibly a loss of derivative in 
the energy). However, here it can be seen that an additional factor of |^| on (I8.17P 
would give a bound of |^| |r/f and thus only two derivatives on the weight, as well 
as a loss of a derivative in L°°. This is what forces us into the splitting into high 
and low frequencies. 

8.4. Adding a Derivative. Once we have established bounds for |ll^x/(t )|^2 
(Propositions 17.21 and 17.3p we can add another derivative: By going back to all 
terms involved we can ascertain that in this argument the bounds for two deriva¬ 
tives on a weight give rise to similar bounds for he. for three derivatives 

on a weight ~ all one needs for this is the product rule of differentiation and the 
bound for the case of two derivatives. In direct analogy to Propositions 17.21 and 17.31 
one then obtains Proposition 17.11 

Its proof consists of revisiting all terms and adding a derivative: By the product 
rule of differentiation this gives terms of the kind we had before with three derivatives 
on the weight instead of two (Regions 1 and 2), terms with two derivatives on the 
weight and one more derivative on the other profile (in the case of Region 3, see 
also Hemark 18.21 above) as well as terms without weights that can be dealt with by 

®The point here is that a possible bound x L°° —>■ just barely fails, so in closing the 
bootstrap for Theorem ITT] we will want to choose p sufficiently small. 
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(higher order) energy estimates. All those can be estimated the way described here, 
the only change being a possible loss of an additional derivative. 
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Appendix A. Proof of the Dispersive Estimate 

Proof of Proposition \5.1\ We localize g in frequency (using Littlewood-Paley the¬ 
ory), estimate the different pieces separately and sum in the end. 

1. Localization: Denoting by Pj the Littlewood-Paley projectors, by Young’s 
convolution inequality we have 


I L\t I 

|e ^1^00 = 






O'GZ 


< J]; \e^^'Pjg\, 

LOO JGZ 


E 






\F-^ [f>{2-^om))\L^ 


ju-i ( 


JGZ 

where Qj is a Littlewood-Paley projector associated to a slightly “fattened” bump 
function (p (which equals 1 on the support of p). 

Thus it suffices to estimate on the frequency localizers p. To this end we 
notice that by a change of variables 

dR2 JR2 

As will be shown further below we have 


(A.l) < Ct ^ 

with a constant depending only on the localizer p. Hence our final estimate will be 


(A.2) 






Ll2-2t 




\Q,g\Lr<^2^^i2-H)-^\Qjg\^, 


j& 




as claimed. 

2. Proof of (jA.ip . We recognize the estimate of 














LONG TIME STABILITY FOR SOLUTIONS OF A /3-PLANE EQUATION 


37 


in L°° as a classical question in the theory of oscillatory integrals and resolve it using 
the method of stationary phase. Here the phase function cf) is given as 

and satisfies 

so that for every (t, x) G M'*' x there are at most four stationary points. A direct 
computation then gives the determinant of the Hessian of 4> as 

4 

det Hm =- 1 

which is non-degenerate away from the origin and bounded above and below in the 
support of ip. Hence the standard stationary phase lemma concludes the proof. □ 

Appendix B. Global Well-posedness (Without Decay) for Large Data 

In this section we prove that the /3-plane equation is actually globally well-posed 
in k > 1, for large data. However, a priori the norms of the solution could 

grow faster than double exponentially in time, which is completely different from 
the situation in the small data case. 

Theorem B.l. For each k > 1 and loq G n H~^, the initial value problem 

\ dtUJ + u-Vuj = Liuj, 

^ ’ \ w(0) = wo G n H-^ 

for the original (i-plane equation has a unique solution uj G C{[0,oo), H^) for all 
time. Recall that here u = V-*-(—A)“^a; and Lioj = 

Remark B.2. Note that the H~^ control on w is merely the conservation of the 
velocity held u. 

Remark B.3. The growth given in our proof is of triple-exponential order. However, 
it is possible to improve the proof to give an upper bound of double-exponential 
type - this is done by carefully checking the constants in the Sobolev imbedding 
theorem and, instead of using L^, estimating through for p close to 2. 

The proof basically adapts the standard technique used to obtain the well-posed- 
ness of the 2d Euler equation to our setting. 

Proof. The main idea is to use the conservation of the vorticity (which has roughly 
the same regularity as Vtt) to prove an a priori bound on the norms of the 
solution. This is done using a borderline Sobolev inequality of the following type: 

I Vul^oo < |w| ^oo log (^ 

for any s > 1, since u = V-'“(—A)“^a;. 
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From this one estimates that 

\{u ■ Vu,u})h‘>\ < |Vtt|^oc \uj\]js 

using, for example, the Kato-Ponce inequality and the fact that u = 

This implies the energy estimate 

dt |w |Hs < |w |LOO log ( Mifl ] . 

\Ml^J 

Now, if we had a priori control on uj in L°°, say 

|w(i)lL-> < Mt), 

we would be able to control the estimates for all time using the Gronwall in¬ 
equality as follows: 


dt \oo\hs < \oj\loo log 






CJ 


L°° 






= \ uj\loo logdwlj^s) \u\jjs - |w|ioo log (|w|ioo) \oo\jjs 
< log |a;| ^ |w| 


since 


Hence 


and so 


\H’> 


, 11 , , 
X ■ log — < a; e 0, oo). 
X e 


dt \uj\hs < { A{t) log \uj\hs + - ) \uj\fjs , 


dtlog\uj\Hs < |^A(t)log , 

which shows that \uj\jjs is bounded in L“^(0,oo). 

A priori global bound on oj in L°°. Our equation reads 

dtOJ + u ■ Vuj = (9i(—A)“^a;. 

Recall that (5i(—A)“^ti;, w)x ,2 = 0 and {u ■ Vu],u})l 2 = 0, thus 

\u{t)\L2 = \u!o\l2 

and similarly 

\uit)\L^ = \uo\l2 , 

just by noting that (V“’“(—A)“^ (tt • Vw) ,u)i 2 = 0. 

From (9i(—A)“^a; = U 2 , the Sobolev embedding theorem and the boundedness of 
Riesz transforms on we deduce that 
|9i(-A)-ic 

< \uo\l^ + \^o\l^ + \^\loo ■ 

We may also write the /3-plane equation as: 

dtuj + u ■ Voj = F, 


LOO ^ ~ I^Il3 + \^\l3 < \'^\m + 
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with 

(B.2) 


l-^l^oo < C'(wo) + 


This gives, using the Gronwall Lemma, a global bound on the L°° norm of u. 
Indeed, define the Lagrangian flow map associated to u by 


A 

dt 

$(0, X 




= X. 


Then 


— (uj o <h) = T o 
dt 


We also have \oj o d>|^oo = |w|^oo and \F o which means that 

|w|loc < IwoIloo + / |-F|^oo is)ds. 

Jo 

Combined with (jB.2p we thus obtain an a-priori bound on uj in L“^(0, oo), say 

|w(t)|ioo < A{t), 

which is all that was left to show. □ 


Appendix C. A Bound on / in 
In this section we prove Corollary 12.41 i.e. we demonstrate that in t 


of Theorem 12.11 we have 




the arguments involved are closely 




< e 2 , provided that 


ier‘^o(o 


le framework 
< e. This 

L°° 


is a consequence of the estimates (12.31) applied to the Duhamel formula (16.ip and 


inked to those in Section El 


Proof. Recall the Duhamel expression for /, 


fit, 6 = /o(0 + [ [ r])f{s, C - r])f{s, r])dr]ds. 

Jo Jr^ 

If at least one of the profiles in this expression is localized to frequencies larger than 
some N > 0 we have the bound 


/ 

Jr‘ 


J - ri)f{s,ri)dtj 




for k as in the proof of Theorem 12.II 

Hence we may assume that all frequencies are bounded by > 0. As for the 
weighted estimate the idea is to use the oscillations of the phase <1> to obtain ad¬ 
ditional decay. To this end we recall that divisions of 

Fourier space akin to those in section 18.31 

Firstly, for small rj we directly obtain the bound 


1 ^ 1 ^ I e^"'^m{^,r])f{C-v)fid)dr] 

<N^ 






Loo 


I/I 


H- 


by estimating one of the profiles in and the other in L°° with derivatives. 
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Secondly, when |r/| > ^ we further split the Duhamel term into two pieces by 
writing 


ler 


JiL 


e' mffdrj = 


e^^^mffdrj 


l>Tn5> l€-2f;l>p 


+ 


[ 

J\n 


= -.A 


>^,IS-2r?|<P 


e^^'^mffdr], 


where p > 0 remains to be chosen later. 

The hrst term of these is then amenable to an integration by parts in p, and the 
only thing left to do is to ascertain that the singularities that arise are controlled. 
We obtain 

1 


A = 


l?-2r?|>p 


- p)f{v)dp 


IS 


|V„^| 


-f 






m 


Vp$| 


fi^-r])f{v)dp =:/ 


1 

is 


IS-2'7|>P 


+ 


-f 


>TSn’ l?-2p|>p 


- p)f{p)dp =: II 
- v)Vr,hp)dp =: III. 




Both terms II and III can be estimated in L°° by 


|^|2// 

L°° ’ 

Infill 

N 

< — 
L°° pS 

\riff 

L2 

Ihl'v,/ 


L2 


since as noted above 
domain of integration. 


V„^> 




= ^CM Np 1 |,e - hP \v? in this 


To control I we compute (akin to the remarks in Section [8.Ijl 

(le+ 1,1-=’) 


m Vr 


iVp$r 


< 

rs_/ 


\vf le-277^^1^ 


and 


< 




le-2r?ll^l 

|Vpm| _ |e||e-r7|2|ry|2 ^ \^ - p(^ 


iVp^l 


|77n^-2r7| lei ^ IC-2771 


We further subdivide the region of integration in I: On the one hand this gives 






'l«-2p|>l 




m 


Vp^l 


/(? - v)fiv)dv 


l/liL*: • 






























LONG TIME STABILITY FOR SOLUTIONS OF A /3-PLANE EQUATION 


41 


On 


the other hand we note that ^ |log(/o)l, so estimating 


1^ — 27/1 ^ in L^, one of the profiles in and the other one in L°° in the region 
where 1 > |^ — 2r]\ > p leaves us with the bound 




< 


1 


L°° S 

Collecting the above estimates shows that 


^iv^iiog(p)i 1/1^3 \em) 


d, 


\em) 


i/ i/ 


L2 


+ \f\l^ + s-^N'^ |log(p)| 1 / 1^2 ler /(O \f\H^ ■ 

^oo 


rl-k I r,2 


Inserting now the bounds for / and revisiting the arguments in the bootstrapping 
scheme in Section [3] yields the claim, once we let p = (recall also that |/|^-i = 

□ 


u 


L2, 


Appendix D. Some Computations 

The following claims are the result of straightforward computations. 

About the Null Forms. Recall that and fh{^,p) := m{^,^ — p). 

Lemma D.l. We have 


Vgm(^,?7) = 

1^1 


(D.l) 




- 2 - 




Tj 

V,,m(C,r/) = - 2-^m{C,y), 

\y\ \y\ 


V,,m(^,?7) = -- 




- 2 tI—^777(^,77). 


Some Phase Derivatives. Direct computation gives the following derivatives. 


o .. yl-yl , i^2 - y2? - {^1 - yif 

^ 771 ^ — . ./I ' 


(D.2) 

(D.3) 

(D.4) 


\^-y\ 




|V„$| = 


\y\" 

\^-y\"\yr 


\^-yr 


iv.^l = 


\y - 2^1 \y\ 
\C-yf\C\" 


dl $ = ° 


'Vl 


'm 


, _ . — 8—F T 6- . 

le-^r 1^1® \y\" 



























42 


TAREK M. ELGINDI, KLAUS WIDMAYER 


(D.5) 




^§2 ^ V 2 ^ 
%^r,2^ = 


^ 2(6 - m) [(^1 - rjif - 3(6 - mf] 

1 ^- 6 ' 

2(6 - m) [- 3(6 - + (6 - mf] 
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